function dx=fun(x)
f1=x(1)^2+x(2)^2-1;
f2=0.75*x(1)^3-x(2)+0.9;
dx=[f1;f2];
end